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ABSTRACT 


This  study  provides  some  insight  to  the  differences  between  the  wave- 
induced  flow  quantities  measured  in  a fixed  and  wave- following  transformed 
coordinate  system.  The  wave-induced  quantities  have  different  flow 
characteristics  depending  on  the  relative  magnitude  of  the  mean  free 
stream  speed  with  respect  to  the  wave  celerity.  In  order  to  detect  the 
instantaneous  critical  layer  mechanism  associated  with  the  modified 
Miles'  theory,  a wind  speed  just  above  the  wave  celerity  should  be  used. 
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1.  INTRODUCTION 


The  momentum  and  energy  transfer  from  the  wind  to  water  waves  on 
the  surface  of  the  ocean  is  a major  geophysical  phenomena.  In  the  past 
two  decades  many  efforts  have  been  made  to  gain  a complete  understanding 
of  the  processes  responsible  for  the  generation  and  the  growth  of  the 
water  waves.  Neumann's  wave  spectrum  concept  was  an  important  advance, 
but  the  first  significant  contributions  to  the  dynamical  wave  theory 
were  contained  in  Phillips  (1957)  and  Miles  (1957) . 

Phillips'  (1957)  resonance  mechanism  describes  the  response  of  the 
water  surface  to  the  direct  action  of  air  turbulent  pressure  fluctuations 
and  is  responsible  for  the  initiation  of  the  water  waves.  The  wave 
growth  rate  associated  with  this  mechanism  is  linear. 

Miles'  (1957)  inviscid,  quas i- laminar , instability  mechanism 
includes  the  feedback  of  air  perturbations  due  to  water  waves  to  the 
waves  themselves.  The  weaknesses  of  this  mechanism  are  the  resulting 
infinite  vorticity  at  the  "critical  height"  where  the  mean  wind  velocity 
is  equal  to  the  phase  velocity  of  the  waves  and  the  requirement  of  a 
free-slip  boundary  condition  at  the  interface.  These  difficulties  were 
overcome  later  by  Miles  (1959)  and  Benjamin  (1959)  by  retaining  viscous 
effects  to  form  a "critical  layer"  near  the  critical  height  and  a viscous 
layer  at  the  interface.  However,  there  are  no  obvious  differences  between 
the  results  from  the  inviscid  model  and  from  the  viscous  model.  Hence, 
for  energy  transfer  from  the  wind  to  the  waves,  the  inviscid  model 
appears  to  remain  valid.  The  physical  key  to  Miles'  theory  is  that 
the  component  of  pressure  perturbations  in  phase  with  the  wave  slope, 
this  compounent  being  the  driving  force  for  the  energy  transfer,  is 


present  due  to  the  existence  of  the  critical  layer.  Therefore,  one 
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concludes  that  the  critical  layer  has  an  important  role  in  the  momentum 
and  energy  transfer  from  the  wind  to  the  waves.  The  growth  rate  of 
water  waves  according  to  this  model  is  exponential. 

The  equations  involved  in  both  Phillips'  and  Miles'  processes 
are  linear,  that  is,  the  components  of  pressure  fluctuations  and  of 
water  waves  with  the  same  frequency  and  wave-number  are  linearly 
related.  In  their  models,  nonlinear  effects  produced  by  the  inter- 
actions among  the  components  of  turbulence  and  wave  perturbations  at 
different  frequencies  and  wave-numbers  are  ignored.  Experimental 
results  both  in  the  field  [Snyder  and  Cox  (1966),  Barnett  and  Wilkerson 
(1967),  and  Dobson  (1971)]  and  in  the  laboratory  [Bole  and  Hsu  (1969) 
and  Stewart  (1970)]  showed  that  the  observed  wave  growth  rates  are 
one  order  of  magnitude  greater  than  those  predicted  by  Miles'  theory. 
There  are  other  experimental  measurements  which  lead  either  directly 
or  indirectly  to  similar  conclusions.  This  inconsistency  between 
predictions  and  observations  suggested  a re-examination  of  the  role 
of  turbulence  which  was  originally  neglected  in  Miles'  formulation. 

To  evaluate  the  turbulent  effect,  one  approach  which  has  been 
used  most  frequently  is  to  construct  a numerical  model  for  the  wave 
perturbation  field  by  retaining  the  wave-induced  turbulent  Reynolds 
stresses  in  the  linearized  wave  perturbation  equations.  Unfortunately, 
this  type  of  approach  has  been  greatly  hindered  by  the  associated 
closure  problem.  Several  closure  models  have  been  proposed  (Long, 

1971;  Davis,  1970,  1972;  Townsend,  1972;  Saeger  and  Reynolds,  1971; 
Norris  and  Reynolds,  1975;  and  Gent  and  Taylor,  1976).  None  of  these 
models  led  to  a satisfactory  prediction  of  the  wave  growth  rate.  All 
the  models  were  based  on  ad  hoc  or  phenomenological  assumptions 
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necessitated  by  the  lack  of  accurate  experimental  data  on  the  Induced 
turbulent  Reynolds  stresses.  According  to  Davis  (1970,  p.  730),  "the 
J problem  of  wave  generation  cannot  be  understood  until  the  wave's 

influence  on  the  turbulent  stresses  is  known."  A direct  measurement  of  the 
wave  perturbation  field  including  the  induced  turbulent  Reynolds  stresses 
seems  to  be  an  expeditious  way  to  reveal  the  wave  perturbed  flow  field 
and  hence  to  achieve  successful  closure  modeling. 

A measurement  of  the  induced  turbulent  Reynolds  stresses  was  first 
attempted  by  Kendall  (1970)  and  subsequently  at  Stanford  by  Yu,  et  al. 

(1973)  and  Chao,  et  al.  (1976).  However,  the  first  conclusive  and 
complete  set  of  data  on  the  induced  turbulent  Reynolds  stresses  was 
presented  by  Hsu,  et  al.  (1977).  One  interesting  feature  shown  in  Hsu, 
et  al.  is  that  an  eddy  viscosity  type  of  closure  model  seems  appropriate 
for  the  induced  turbulent  Reynolds  stresses. 

In  addition  to  searching  for  a closure  model,  a complementary 
, approach  is  to  evaluate  directly  the  turbulent  effect  on  the  Miles' 

mechanism.  When  turbulence  exists,  the  role  of  viscosity  in  smoothing  the 
concentrated  vorticity  distribution  in  the  critical  layer  is  taken  over 
by  the  more  diffusive  turbulent  mixing.  The  physical  significance  of 
the  turbulent  effect  on  the  critical  layer  structure  was  first  acknowledged 
by  Lighthill  (1962)  and  then  followed  up  by  Stewart  (1967,  1974). 

According  to  Lighthill,  the  turbulence  effect  was  supposed  to  broaden  the 
critical  layer  and  to  render  Miles'  mechanism  more  efficient.  But,  the 
extent  of  the  enhancement  of  the  energy  transfer  from  the  wind  to  the 
waves  due  to  the  turbulence  was  not  estimated  by  Lighthill.  The 
Hasselmann's  (1968)  theory  also  suggested  that  the  turbulence  will  cause 
the  critical  layer  to  extend  over  a wide  range  of  the  wind  velocity  profile. 
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The  turbulent  effect  on  the  wave  generation  processes  is  twofold: 

(a)  to  transfer  directly  the  energy  from  wind  to  waves  through  the  induced 
turbulent  Reynolds  stresses,  and  (b)  to  enhance  the  Miles'  mechanism  by 
changing  the  characteristics  of  the  critical  layer.  Because  the  energy 
contained  in  the  turbulence  is  relatively  small  compared  with  the  amount 
of  energy  required  for  the  wave  growth,  Stewart  (1967)  inferred  that 
transfer  (a)  should  be  insignificant.  He  also  interpreted  transfer  (b) 
in  terms  of  the  so-called  "catalytic  effect".  The  concept  of  the  catalytic 
effect  was  further  extended  by  Stewart  (1974)  through  examining  and 
comparing  the  transfer  by  Miles ' mechanism  while  using  either  the  mean 
velocity  profile  or  the  instantaneous  velocity  profiles,  as  measured  by 
Kim,  et  al.  (1971).  The  visual  studies  of  Kim,  et  al.  suggested  that 
approximately  90  percent  of  the  turbulent  Reynolds  stress  is  produced 
during  a wave-like  "bursting"  process.  Their  measured  turbulent  Reynolds 
stress  was  also  consistent  with  their  prediction  based  on  the  measured 
instantaneous  velocity  profiles.  Hence,  it  seems  likely  that  the  wave- 
associated  Reynolds  stress  can  be  regarded  as  the  statistical  average  of 
those  stresses  produced  by  the  critical  layer  mechanism  according  to  the 
instantaneous  velocity  profiles.  These  concepts  and  the  philosophy 
of  Davis  (1974,  p.  673),  that  "all  averaging  is  delayed  until  the  dynam- 
ical equations  have  been  solved  ...",  has  enabled  us  to  consolidate  the 
turbulent  catalytic  effect  into  mathematical  expressions.  The  result  is 
then  regarded  as  a modified  Miles'  theory. 

The  experimental  aspect  of  this  research  program  contains  two  facets. 

One  is  the  study  of  the  structure  of  the  wave-induced  flowfield,  expecially 
for  the  wave-induced  turbulent  Reynolds  stresses  which  form  the  kernel 
of  the  closure  problem.  The  other  is  the  study  of  the  turbulent  effect  on 
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the  structure  of  critical  layer  as  it  appears  in  our  modified  Miles' 
theory;  "he  details  of  which  are  given  in  Section  2.  For  the  first 
facet,  we  extended  the  measurements  of  Rsu,  et  al.  (1977),  which  were 
performed  at  a mean  free  stream  velocity  of  2.4  m/sec  and  in  a wave- 
following, transformed  coordinate  system,  to  mean  free  stream  velcolties 
of  1.4  and  2.9  m/sec  and  to  the  inclusion  of  a fixed  frame  measurement 
to  examine  the  difference  between  fixed  frame  and  transformed  frame 
measurements  under  identical  experimental  conditions.  For  the  second 
facet,  we  measured  the  instantaneous  velocity  gradient  by  using  two  X-array 
hot-film  probes  separated  vertically  by  1.2  cm.  The  statistical  properties 
of  the  instantaneous  velocity  gradient  are  the  key  to  the  evaluation  of 
the  modified  Miles'  theory  as  shown  in  Section  2.  The  facility  used 
for  the  experiment  is  described  in  Section  3;  the  data  acquisition  and 
reduction  are  then  given  in  Section  4.  In  Section  5 we  describe  the 
experimental  results,  followed  by  a concluding  remark  (Section  6). 


2.  THEORETICAL  ASPECTS 


2.1  Momentum  and  Energy  Transfer  from  Wind  to  Waves  — a Modified  Miles' 
Theory 

Consider  a turbulent  wind  flowing  over  a sinusoidal  water  wave  which 
propagates  in  the  x-direction  and  is  given  by  rf  = a • cos(kx  - u)t)  , 
where  a is  the  wave  amplitude,  k is  the  wave-number  and  u)  is  the 
circular  wave-frequency.  The  vorticity  equation  in  the  lateral  (z-direction) 
can  be  obtained  by  eliminating  the  pressure  terms  in  the  momentum  equations 
in  the  x-  and  the  (vertical)  y-directions.  The  resultant  equation  is 


where  V is  the  kinematic  viscosity,  = (3^/9x^‘)  + (3^/9y^)  + (9^/3z^) 
and  u , v , and  w are  the  velocity  components  in  the  x , y , and  z 
directions  respectively.  In  (2.1)  we  have  assumed  that  the  air  is 
imcompressible  so  that  the  continuity  equation  is 


9u  9v  9w 
3x  9y  9z 


0 


(2.2) 


We  assume  that  the  mean  flow  is  parallel,  i.e.,  u = u(y)  , 
v = 0 , and  w = 0 , and  that  the  wave-induced  flow  is  two  dimensional 
in  (x,  y) . Hence,  we  have 

u = u(y)  + u(x,y,t)  + u'(x,y,z,t) 

v = v(x,y , t)  + v’(x,y,z,t)  (2.3a,b,c) 

w = w' (x,y,z,t) 

-6- 
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Here  denotes  a wave-induced  component  and  " ' " denotes  turbulent 
fluctuations,  while,  in  order  to  separate  the  mean  and  the  wave-induced 
quantities  from  the  background  turbulence,  we  define  a time  average  as 


i f m 

g = lim  - J g(x,  t)  dt 
T -*■<*>  -T/2 


(2.4) 


and  a phase  average  as 


<8>  = lim  2nVt  S 8(-’  t+nx) 

N-*-°°  n = -N 


(2.5) 


where  g is  a flow  quantity  of  interest,  x is  the  position  vector  at 
which  the  averages  are  performed  and  x is  the  period  of  wave  oscillation. 
The  wave-induced  flow  then  can  be  found  from 


g = <8>  ~ g 

The  vorticity  in  the  z-direction  is  given  by 


n = ly  " lx  = + y*  ^ + y»  z» 


(2.6) 


(2.7) 


Hence, 


o’  _ _9u  n - n'  - 3u'  3v' 

3y  * 9y  3x  ’ “ “ 3y  " 9x 


(2.8a,b,c) 


by  substituting  (2.3a,b,c)  into  (2.7). 

When  (2.3a,b,c)  and  (2.7)  are  substituted  into  (2.1),  the  resultant 
equation  is  first  time  and  then  phase  averaged  and  the  time-averaged 
result  is  subtracted  from  the  phase-averaged  result,  we  find 
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It  has  been  shown  by  Davis  (1969)  and  Hsu,  et  al.  (1978)  that  s^  is 

related  to  the  nonlinear  effect  of  the  wave-induced  motion  and  irrelevant 

to  the  momentum  and  energy  transfer  from  wind  to  waves,  although  the  wave- 

induced  flow  structure  may  be  different  for  a large  amplitude  wave. 

Under  the  small  amplitude  assumption  on  the  water  wave,  s^  can  be 

2 

safely  neglected,  since  it  is  only  of  order  (ka)  . On  the  other  hand, 
the  terms  with  r\^  in  (2.9)  are  of  order  (ka)  and  are  assumed  to  be 
significant  in  changing  the  wave-induced  flow  structure  relevant  to  the 
momentum  and  energy  transfer  from  wind  to  waves.  However,  the  inclusion 
of  r\j  in  (2.9)  raises  the  well  known  "closure"  problem  because  the 

are  the  additional  unknowns  which  make  the  problem  unsolvable  unless 
some  closure  relations  are  assumed.  The  viscous  terms  in  (2.9)  are 
usually  negligible  when  r^  is  included  in  (2.9),  except  near  the 
interface  where  the  turbulence  is  damped  out  by  viscosity. 
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It  must  be  pointed  out  that  the  application  of  the  averaging 
procedures  here  has  simplified  the  problem  by  artifically  restricting 
the  critical  layer  in  a mean  sense  to  a region  near  the  definite  position 
where  u - c = 0 . Furthermore,  it  also  forced  us  to  rest  the  solution 
of  (2.9)  entirely  on  the  closure  relations  associated  with  the  wave- 
induced  turbulent  Reynolds  stresses  r*  . A successful  closure  model 
should  not  only  lead  to  a description  of  the  effect  of  turbulence  in 
changing  the  structure  of  critical  layer  as  related  to  the  wave-induced 
motion,  but  also  to  a prediction  of  precisely  the  amount  of  momentum  and 
energy  transfer  from  wind  to  waves.  Unfortunately,  to  date,  there 
exists  no  such  applicable  model. 

It  is  important  to  note  that  in  a real  flow  the  effects  of  non- 
linearity and  of  turbulence  will  change  the  location  of  the  critical  layer 
from  time  to  time  according  to  u + u'+u'  - c = 0.  Thus,  the  slope  and 
the  curvature  at  the  critical  height  at  any  time  are  strongly  influenced 
by  these  effects.  Stewart  (1974)  suggested  that  the  mean  velocity 
profile  was  not  the  appropriate  thing  and  that  consideration  of  the 
instantaneous  velocity  profile  would  be  helpful  in  describing  the  energy 
transfer  from  wind  to  waves.  In  order  to  keep  the  spirit  of  his  "catalyst 
effect,"  after  substituting  (2.3a,b,c)  into  (2.1)  and  subtracting  the  time 
average  of  the  resultant  equation  from  the  resultant  equation,  the 
difference  can  be  written  as 


(u  - c)  -jp  - T’  + f - 
dx  n 


v^±+  f + f 

a 2 v 

3y 


(2.12) 


. aa  an 
When  9F  = " C 9^ 


organized  component 
etc.:  and  f 


3u 
u -r— 

3y 


is  applied.  In  (2.12),  T’  is  the  term  containing  no 
, i.e.,  <T’>  =0;  f^  is  the  nonlinear  term  containing 

and  fv  are  expressed  as 


^9x2  3y2 


\ ~ 3 /' 

^ r12  3x3y 


11  22 


(2.13) 


r - v (4  * 4\ 

V W 3y2/ 


(2.14) 


respectively. 


Since,  we  are  interested  in  the  role  of  the  instantaneous  velocity 

profile  in  the  framework  of  traditional  Miles'  theory,  T'  and  f in 

n 

»v 

(2.12)  can  be  ignored.  The  solution  for  in  the  limit  of  very  high 
Reynolds  number  can  be  found  from  (2.12)  by  letting  f^  be  zero  while 
admitting  a singularity  at  y ^ where  u - c = 0 as  generally  imposed 
in  the  inviscid  theory.  To  the  first  order  approximation  in  ka  , fi  , 
v , and  f are  sinusoidal  functions  of  x . By  dividing  (2.12)  by 
(u  - c)  and  taking  the  phase  average,  we  have 


WHx,  y.  t)  - H (x  - £ , y,  t) 


+ i/(u-  c)  * f ~ y,  t) 


(2.15) 


Note  that  in  reaching  (2.15),  the  dynamical  equation  for  the  vorticity  is 
solved  before  the  phase  average  is  taken,  as  emphasized  by  Davis  (1974). 

Hsu,  et  al.  (1977)  showed  that  the  momentum  and  energy  transfer 
from  wind  to  waves  is  mainly  produced  by  the  wave-induced  pressure  which 
is  supported  by  the  wave  form  to  give  momentum  flux  and  acts  against  the 
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vertical  oscillation  to  do  work.  The  momentum  and  energy  transfer  by 

the  induced  pressure  can  also  be  interpreted  in  terms  of  the  wave- 

associated  Reynolds  stress.  According  to  Lighthill  (1962),  the  wave- 

associated  Reynolds  stress  T (y)  is  produced  from  the  transport  of  the 

w 

wave  induced  vorticity  f)(x,y,t)  by  the  wave-induced  vertical  velocity 
v(x,y,t)  , i.e. , 


Tw(y)  = - p uv  (y)  = p J ftv  (£)  d£  (2.16) 

y 

Because  only  the  component  of  fi(x,y,t)  in  phase  with  v(x,y,t) 
can  correlate  with  v to  produce  the  wave-associated  Reynolds  stress, 
the  first  term  on  the  right  hand  side  of  Eq.  (2.15)  contributes  only  in 
the  critical  layer.  However,  the  second  term  can  contribute  either  when 
u = 0 or  u r 0 depending  on  the  phase  of  f relative  to  v . 

Substituting  Eq.  (2.15)  into  Eq.  (2.16)  and  applying  the  residual  theorem 
of  complex  analysis,  we  find 


T (y)  = p 

w 


+ p 


|l  - ct(S)]  ‘ (2.17) 


where  a(£)  is  the  probability  density  of  the  occurrence  of  the  critical 

layer  in  the  interval  | d^  , £ + d?)  • Wien  Eq.  (2.17)  is 

evaluated  at  y = 0 (actually  outside  the  viscous  sublayer  because  the 

viscosity  is  ignored),  x (0)  should  give  the  momentum  transferred  to 

w 


2.2  Transformed  Wave-Following  Coordinate  System 


It  was  shown  by  Hsu,  et  al.  (1977)  that  the  mean  flow  follows  the 


wave  form.  Hence,  when  the  critical  layer  is  low  as  usually  encountered 


in  a wave  growth  period,  most  of  the  instantaneous  velocity  profiles  have 


their  critical  layer  below  the  wave  crests.  It  is  then  difficult  to 


collect  all  the  information  in  a fixed  frame.  The  difficulty  of  using  a 


fixed  frame  can  be  overcome  by  using  a transformed  wave-follwoing  frame 


given  by 


t*  = t 


x*  = x 


(2.18a,b,c) 


y*  = y - f n 


where  n is  the  displacement  of  the  wave  from  mean  water  level. 


f(v*)  = ■_sinh.(kH  - ky*) 
sinh(kH) 


(2.19) 


and  H is  the  height  from  the  mean  water  level  to  the  roof  of  a wind. 


wave  channel.  Thus  y*  = 0 corresponds  to  the  interface  y = h , and 


y*  = H corresponds  to  y = H , the  channel  roof.  The  lines  of  constant 


y*  are  the  streamlines  of  a potential  wave  motion  in  the  absence  of  wind. 


The  adoption  of  the  transformed  coordinate  system  allows  a complete 


description  of  the  flowfield.  The  difficulty  of  locating  and  dealing 


with  an  undulating  critical  layer  over  the  wave  form  is  also  nicely 


circumvented. 


In  the  transformed  coordinate  system,  the  wave-induced  velocities 


u(x*  , y*)  and  v(x*  , y*)  are  regarded  as  those  organized  oscillatory 


velocities  around  u(y*)  , which  is  the  mean  velocity  along  a curve 


■ 


rT 


where  y*  = constant  . The  substantial  difference  between  7T(x*  , y*) 
and  “u(x  , y)  was  shown  by  Hsu,  et  al.  (1977).  They  also  showed  that 
u*(x*  , y*)  and  v*(x*  , y*)  defined  by 


u*(x*  , y*)  = u(x*  , y*)  - fn 


v*(x*  , y*)  a v(x*  , y*) 


(2. 20a,b) 


have  an  equivalent  role  to  u(x,y)  and  v(x,y). 

The  difference  in  the  vorticity  in  the  two  coordinate  systems  is 
noteworthy  since  the  wave-induced  vorticity  is  important  in  the  pro- 
duction of  the  wave-associated  Reynolds  stress.  The  vorticity  of  the 
flow  is  given  by 


3v 

3x 


3u 

3y* 


3f  ~ 3u  3v 
3y*  ^ 9y*  " 9X*  + 


f 


3n  3v 

3x  3y* 


+ o(kn)2 


(2.21) 


Taking  a time  average  (or  average  along  x*)  , we  have  to  the  first  order 
in  krf 


n<y*) 


3u 

3y* 


(2.22) 


Hence,  the  mean  vorticity  in  the  two  frames  remains  the  same  because 
the  mean  velocity  u in  the  two  frames  is  the  same  up  to  the  first 
order  in  kn  (Hsu,  et  al.,  1977).  The  wave-induced  vorticity  fi(x*,y*) 
measured  in  (x*,y*)  is  then  obtained  by  subtracting  Eq.  (2.22)  from 
the  phase  average  of  Eq.  (2.21)  and  is 


fl(x*,y*) 


3vT  3f  ~ 3u  _ 3v 

3y*  3y*  3y*  3x* 


(2.23) 


.1 

3 
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The  vorticity  fl*(x*,y*)  is  defined  by 


ft*(x*,y*) 


When  Eqs.  (2.20a,b)  and  (2.23)  are  introduced  into  Eq.  (2.24),  we  find 


(J*(x*,y*)  = fi(x*,y*)  - frT 


which  has  a form  similar  to  Eq.  (2.20a) 


The  wave-associated  Reynolds  stress  in  the  transformed  coordinate 


system  is  now  given  as 


T (y*)  = -p  “u*v*  (y*)  = p 


We  also  have 


u*v*  (y*)  = u v(y*)  - f 


3.  EXPERIMENTAL  APPARATUS  AND  INSTRUMENTATION 


3.1  The  Air-Water  Channel 

The  channel  used  for  the  experiments  was  described  in  detail  by 
Hsu  (1965)  and  subsequent  investigators  at  Stanford.  It  was  designed 
to  facilitate  the  generation  of  water  waves  by  wind  and/or  by  a mechanical 
wave  generator  for  the  experimental  study  of  the  wave  generation  problem. 
Figure  1 shows  a schematic  of  the  channel.  It  is  38  m long,  1 m wide  and 
1.9  m high.  In  this  study,  the  distance  between  the  channel  roof  and  the 
mean  water  level  is  0.97  m.  The  glass-walled  test  section  is  approxi- 
mately 24  m long,  remaining  from  the  air  entrance  to  the  down-stream 
beach.  The  data  taking  station  was  located  13  m from  the  air  entrance. 

Wind  is  produced  by  drawing  air  through  the  test  section  with  a 
suction  fan  at  the  downstream  end  of  the  channel.  The  wind  field  was 
conditioned  by  a thick  honeycomb  and  a group  of  fine-meshed  screens 
installed  at  the  exit  of  the  air  inlet  and  by  a second  honeycomb 
located  just  in  front  of  the  fan.  To  extend  the  experiment  performed 
by  Hsu,  et  al.  (1977)  which  was  made  at  a mean  free  stream  velocity,  U^,  of 
2.4  m/sec,  the  mean  free  stream  velocities  selected  in  this  study  were 
1.45  and  2.92  m/sec. 

The  water  wave  generator  is  a horizontal-displacement  oscillating 
plate  driven  by  a hydraulic  cylinder  and  subjected  to  a closed-loop 
servo-control  system.  The  wave  plate  is  situated  well  upstream  from 
the  entrance  of  the  test  section  to  form  a forebay,  which  acts  as  a 
transition  region  for  the  development  of  mechanically-generated  water 
waves.  A sloping  beach,  formed  by  a series  of  rectangular  baskets 
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filled  with  stainless  steel  turnings,  is  located  at  the  downstream  end 
of  the  channel  to  reduce  the  wave  reflection.  The  reflected  waves 
amplitudes  were  less  than  5 percent  of  the  incident  mechanically- 

t 

generated  wave  amplitudes.  The  frequency  of  the  water  wave  used  in  this 
study  was  1 Hz.  The  wave  length  determined  from  the  dispersion  relation 
of  deep  water  waves  was  1.56  m and  the  wave  number  was  4.025  m 
Hence,  the  wave  speed  is  1.56  m/sec  which  is  between  the  two  selected 
mean  free  stream  velocities. 

3.2  Wave-Follower  System 

The  wave-follower  system  developed  by  Yu,  et  al.  (1971)  and  modified 
by  Hsu,  et  al.  (1977)  was  used  as  the  primary  instrument  for  the 
measurements.  It  consists  of  mechanical  and  electrical  systems  and  a 
wave-follower  elevator.  The  mechanical  portion  contains  a low- inertia 
motor  and  a vertical  motion  mechanism  that  has  an  aluminum- channel 
holding  a stainless-steel  tube  guided  by  nylon  bushings.  The  tube  is 
driven  by  the  motor  through  a pulley-and-cable  assembly.  A Pitot- 
static  tube  and  two  x-array  hot-film  probes  separated  vertically  by 

1.2  cm  were  mounted  on  a plexiglass  frame  that  was  attached  to  the 
stainless-steel  tube.  The  electrical  portion  is  a control  panel  that 
implements  the  wave-follower  by  controlling  the  motion  of  the  stainless- 
steel  tube.  The  wave-follower  system  is  mainly  controlled  by  two  dials 
on  the  control  panel:  one  pre-sets  the  distance  from  the  mean  water 
level  to  the  probes  and  the  other  determines  the  amplitude  of  the 
oscillation  of  the  stainless-steel  tube  and  the  attached  sensors.  The 
first  dial  gives  the  preset  distance  accurate  to  +0.25  mm.  The  second 


F 
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dial  setting  depends  on  the  calibration  of  a wave  height  gauge  whose 
output  signals  was  used  to  guide  the  wave  follower.  The  wave-follower 
elevator  extends  the  traverse  range  of  the  wave-follower  to  about  60  cm 
to  cover  the  whole  boundary  layer.  The  elevator  was  seated  on  top  of 
the  channel  roof.  The  wave-follower  was  mounted  on  a moving  plate  of 
the  elevator  with  the  aluminum  channel  and  the  stainless  steel  tube 
extending  into  the  wind,  wave  channel.  Figure  2 shows  the  arrangement 
of  the  components. 

3.3  Instrumentation  and  Calibration 

The  sensors  consist  of  a wave  height  gauge  mounted  on  an  aluminum- 
angle  support  fixed  to  the  wind,  wave  channel  roof  and  the  Pitot-static 
tube  and  the  two  x-array  hot-film  probes  attached  to  the  lower  end  of 
the  wave  follower.  The  arrangement  of  the  measuring  probes  is  shown  in 
Figure  2.  Because  the  wave  height  sign5*!  was  used  as  input  for  the 
wave-follower  system,  the  sensing  wire  of  the  wave  height  gauge  was 
installed  approximately  4 cm  upstream  of  the  other  sensors  to  correct 
for  the  phase  lag  of  the  wave-follower  system.  The  wave  height  gauge, 
the  Pitot-static  tube,  and  the  two  hot-film  probes  were  laterally 
separated  by  approximately  5 cm  to  avoid  interference. 

The  wave  height  gauge  measured  the  elevation  of  the  instantaneous 
water  surface.  The  gauge  is  a capacitance  type  and  is  made  of  No.  32 
Nylclad- insulated  coper  wire  stretched  between  the  arms  of  a U-shaped 
frame.  A complete  description  of  the  gauge  can  be  found  in  Colonell 
(1966).  The  calibration  of  the  wave  gauge  was  carried  out  by  mounting 
the  gauge  on  a vertical  traversing  mechanism  which  varied  the  depth  of 
submergence  of  the  wire  in  still  water.  The  accuracy  of  the  calibration 
was  +0.25  mm.  Figure  3 shows  the  calibration  curve. 
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The  Pitot-static  tube  was  used  to  cross  check  the  mean  velocity 
obtained  from  the  hot-film  probes  and  to  calibrate  the  hot-film  probes. 


The  outer  diameter  of  the  tube  is  2.34  mm.  Two  leads  for  the  total 
pressure  and  the  static  pressure  from  the  tube  were  connected  to  a Pace 
differential  transducer  (Model  P90D).  The  Pace  transducer  was  calibrated 
against  a Combust  micromanometer,  with  resolution  of  0.006  mm  of  silicone 
oil  (specific  gravity  0.82).  Preliminary  mean  velocity  data  showed  the 
Pitot-static  tube  and  the  hot-film  probe  results  agreed  to  within  3 per- 
cent. Hence,  the  signal  from  the  Pitot-static  tube  was  not  recorded 
in  the  final  data  taking. 

The  use  of  two  X-array  hot-film  probes  was  to  measure  the 
instantaneous  values  of  3u/3y  and  to  obtain  the  wave-induced  flow 
components.  The  two  hot-film  probes  were  separated  vertically  by  1.2 
cm.  Two  different  probes  were  used.  The  upper  probe  is  the  (old) 
probe  used  by  Hsu,  et  al.  (1977)  and  the  lower  probe  is  a mini-probe 
(TSI  Model  1276) . The  old  probe  was  connected  to  a TSI  1010A  constant 
temperature  anemometer  (also  used  by  Hsu,  et  al.)  and  the  new  probe 
to  one  port  of  a set  of  newly-built  anemometers**.  This  arrangement 
allows  us  to  cross  check  the  performance  of  the  new  hot-film  anemometer 
system.  Calibration  of  the  hot-film  probes  was  carried  out  by  locating 
the  probes  in  the  core  region  of  air  flow  and  changing  the  wind  speeds. 


The  original  plan  was  to  measure  the  instantaneous  velocity  profile  by 
using  ten  miniature  probes.  Hence,  the  newly-built  anemometers 
consist  of  20  channels  of  electronic  circuitry  for  ten  crossed  hot-film 
probes  and  the  necessary  power  supplies.  Each  channel  has  a constant 
temperature  anemometer,  a four-pole  butterworth-type  low  pass  filter,  a 
zero  offset,  and  an  active  amplifier.  Unfortunately,  the  measurement 
of  the  instantaneous  velocity  profiles  was  not  completed  due  to  the 
termination  of  this  project. 
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using  the  Pitot-static  tube  velocity  as  a reference.  Figure  4 shows  the 
comparison  of  the  typical  calibration  curves  between  the  old  and  the  new 
systems.  It  is  seen  that  the  new  system  has  almost  the  same  sensitivity 
(characterized  by  the  slopes  of  the  curves)  as  the  old  system;  however, 
the  new  system  has  a lower  output  than  the  old  system  as  seen  from  the 
output  power  of  the  two  probes.  The  hot-film  probes  were  calibrated 
in  situ  immediately  prior  to  data  taking.  The  uncertainty  in  the  hot- 
fi.lm  probe  calibration  was  mainly  caused  by  the  uncertainty  in  the  Pitot- 
static  tube  and,  hence,  was  approximately  3 percent. 
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4.  DATA  ACQUISITION  AND  REDUCTION 

The  wave  height  and  velocity  data  were  taken  simultaneously  and 
recorded  by  a data  acquisition-reduction  system.  The  system  consists 
mainly  of  an  HP  2100A  computer,  an  HP  2313  analog-to-digital  converter, 
an  HP  7970B  nine-track  magnetic  tape  drive  and  peripheral  input-output 
devices  (Takeuchi  and  Mogel,  1975).  All  the  signals  were  zero- 
suppressed,  amplified  and  low-pass  filtered  at  500  Hz  to  fulfill  the 
Nyquist  criterion  as  samples  were  taken  every  0.001  seconds  for  three 
minutes.  The  sampled  data  were  stored  on  digital  magnetic  tape. 

In  summary,  we  have  two  sets  of  data  at  mean  free  stream  wind 
velocities  of  1.4  and  2.9  m/sec.  The  mechanically-generated  water  wave 
frequency  is  1 Hz . The  amplitudes  of  the  mechanically-generated  water 
waves  for  the  1.4  and  2.9  m/sec  wind  speeds  were  2.88  and  2.67  cm, 
respectively.  Each  run  consists  of  velocity  data  taken  at  30  elevations 
ranging  from  2.0  cm  to  41.0  cm  above  the  interface  plus  the  wave  height 
data. 

The  time  average  and  the  phase  average  defined  by  (2.4)  and  (2.5) 
are  the  two  main  schemes  used  for  data  reduction.  The  wave-induced 
components  were  then  deduced  from  (2.6).  To  obtain  the  amplitude  and 
the  phase  of  the  wave- induced  components,  the  wave- induced  components 
deduced  from  (2.6)  were  analyzed  by  a fast-Fourier-transform  routine 
to  obtain  the  cross-spectral  density  between  the  wave-induced  component 
and  the  surface  water  wave  rf . The  details  of  the  fast-Fourier- 
transform  routine  can  be  found  in  Bendat  and  Piersol  (1971,  Chapter  9). 
The  detailed  data  reduction  scheme  was  given  in  Hsu,  et  al.  (1977). 
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In  deducing  the  statistical  behavior  of  the  instantaneous  velocity 


gradient,  the  instantaneous  velocities  of  the  two  probes  were  reduced 
first  and  then  the  occurrance  of  the  critical  layer  between  the  two 
probes  was  determined  by  finding  if  the  velocity  of  the  upper  probe  is 
greater  than  the  wave  celerity  and  if  the  velocity  of  the  lower  probe 
is  less  than  the  wave  cererity,  or  vice  versa.  The  probability  density 
a(y)  is  then  obtained  by  dividing  the  number  of  the  occurrances  of  the 
critical  layer  between  the  two  probes  by  the  total  number  of  the  data 
examined . 


5.  EXPERIMENTAL  RESULTS 


The  reduced  data  consist  of  those  for  the  mean  air  flow,  the  wave- 

induced  flow  and  the  probability  distribution  a(y)  . The  flow  quantities 

are  usually  presented  in  profile  distributions  as  a function  of  y . 

Because  for  a fixed  frame,  y is  a constant  and  is  equal  to  the  average 

elevation  y*  of  a transformed  wave-following  frame  at  the  same  data 

taking  station,  y is  used  in  abscissa  of  the  figures  to  be  presented 

below.  The  profiles  are  usually  in  non-dimensional  form;  U and 
2 

Uro  are  used  to  normalize  the  velocities  and  the  Reynolds  stresses, 
respectively.  The  abscissa  y is  then  normalized  by  1/k  . The  wave- 
induced  quantity  “g  is  generally  expressed  as 


i (kx-wt ) 
e 


+ conjugate]  + harmonics 


= I g I • cos  (kx  - wt  + 0~)  + harmonics 

8 


where  |g|  is  the  amplitude  and  0~  is  the  phase  lag  angle  of  the 

6 

fundamental  mode.  The  phase  lag  is  with  respect  to  time  using  ri  as 
a reference.  In  this  study,  the  wave-induced  quantity  can  be  approximated 
by  its  fundamental  mode  without  changing  the  overall  conclusions,  because 
the  harmonics  are  relatively  weak. 


; 


frames  is  not  noticably  different.  This  is  the  case  because  as  proven  by 


Hsu,  et  al.  (1977)  their  difference  should  be  to  the  second  order  in  ka 


A wake  characteristic  near  the  free  stream,  as  also  observed  by  Hsu,  et  al. 


(1977),  is  evident  in  the  profiles.  The  friction  velocities  determined 


from  the  profile  method  are  5.36  and  10.6  cm/sec  for  U = 1.45  and  2.92 


2 2 -4 

m/sec  respectively.  Hence,  the  drag  coefficient  u^  /U^  * 13.0  • 10 


which  agrees  with  the  value  obtained  by  Hsu,  et  al.  (1977),  Figures  7 


and  8 show  the  distributions  of  the  mean  turbulent  Reynolds  stresses, 


Apparently,  u'u'  is  one  order  in  magnitude  greater  than  -u'v'  and 


2 2 2 

v'v'  . Reasonable  agreement  between  -u'v'/U  and  u.  /U  is  found. 

00  X 00 


The  lower  u^uj  for  the  probe  B than  that  for  probe  A in  Figure  7 may 


be  due  to  the  misalignment  of  probe  B to  the  flow  direction.  Because  the 


hot-film  sensor  in  probe  B was  broken  after  finishing  the  first  run 


(U^/c  * 0.93).  The  second  run  (U^/c  = 1.87)  was  performed  with  a new 


hot-f,.lm  sensor  for  probe  B.  The  new  probe  B has  a better  alignment 


with  probe  A;  hence,  they  produce  little  difference  as  shown  in  Figure  8. 


5.2  Wave-Induced  Flowfield 


For  the  wave  induced  flowfield,  we  will  present  the  data  for  the 


wind  velocities  at  1.45  and  2.92  m/sec  separately,  because  these  two  sets 


of  data  belong  to  different  flow  regimes,  namely,  U /c  = 0.93  < 1 and 


U /c  = 1.87  > 1 . We  will  see  that  the  characteristics  of  the  wave- 


induced  flow  for  these  two  flow  regimes  are  quite  different.  This  may  be 


because  there  is  no  critical  layer  in  the  first  flow  while  with  higher 


wind  speed  a critical  layer  does  exist.  The  importance  of  the  critical 


layer  is  then  revealed.  We  also  find  that  the  momentum  and  energy  transfer 


in  the  first  case  is  from  waves  to  wind  while  in  the  second  case  is  from 


wind  to  waves. 


* 


5.2a  No  Critical  Layer  Case  (U^/c  = 0.93  < 1) 

Figures  9 and  10  show  the  amplitude  and  phase  for  the  wave-induced 
velocities  u and  v respectively.  The  difference  between  the  results 
obtained  in  the  fixed  frame  and  the  transformed  wave-following  frame 
is  obvious.  While  the  amplitude  |u|  measured  in  the  fixed  frame  is 
considerably  greater  than  that  measured  in  the  wave-following  frame, 
the  amplitude  of  v in  the  fixed  frame  results  is  considerably  lower 
than  the  transformed  frame  results.  The  higher  iT  in  the  fixed  frame 
is  apparently  due  to  the  behavior  of  the  mean  flow  which  follows  the 

wave  form  so  that  the  fixed  probe  picks  up  a component  fTi3u/9y  . 

However,  we  find  no  explanation  for  the  difference  between  tha  two 

frame  results  in  the  amplitude  [v|  ; this  phenomena  is  also 

observed  by  Chao,  et  al.  (1976).  Further  investigations  on  this  aspect 

is  required.  The  phase  0~  is  approximately  180°  as  expected  from  the 

potential  flow  analysis.  The  values  of  0~  increase  to  210°  as  the 

interface  is  approached  but  have  a tendency  to  return  back  to  180°  near 

the  interface.  The  phase  of  0~  is  almost  constant  at  270°  as  it  should 

v 

be  according  to  the  interface  boundary  condition.  In  general,  the 
wave-induced  flow  at  this  wind  speed  is  similar  to  the  flow  inferred 
from  a potential  flow  analysis,  with  minor  modification  on  the  amplitude 
distribution  by  turbulence.  Again,  the  deviation  between  the  results 
obtained  by  probes  A and  B in  Figure  10(a)  may  be  due  to  the  misalignment 
of  probe  B. 

Figure  11  shows  the  wave-associated  Reynolds  stress  -uiv  versus 

' 

ky  . Near  the  interface  (but  outside  a viscous  Stokes  layer),  the  wave- 
associated  Reynolds  stress  is  negative.  This  means  that  the  momentum 
and  energy  transfer  is  from  the  waves  to  the  wind.  That  is  that  the 
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wind  is  lagging  behind  the  propagation  of  the  waves  so  as  to  drag  on 


the  waves  and  result  in  their  attenuation.  In  Section  2.2,  we  have 

shown  that  the  values  of  -*uv  measured  in  the  two  frames  should  differ 

by  f3u/3y*riv.  Since  rf  and  v are  almost  90°  out  of  phase,  n v 

is  probably  very  small.  As  a result,  no  significant  difference  in  -uv 

between  the  two  frames  is  found  in  Figure  11. 

Figures  12,  13  and  14  show  the  amplitudes  and  phases  for  the  wave- 

induced  turbulent  Reynolds  stresses  as  functions  of  ky  . Hsu,  et  al. 

(1978)  showed  that  is  approximately  twice  r^  and  is  one  order 

of  magnitude  greater  than  r^  near  the  interface.  This  is  the 

situation  for  r^  , r ^ and  r*^  shown  in  Figures  12,  13  and  14. 

In  these  figures  we  also  find  no  phase  jump  associated  with  the  phase 

distributions  of  0~^  , ®r^2  an<^  ^22  " -*-s  different  from  what 

was  observed  by  Hsu,  et  al.  (1977)  for  U^/c  = 1.57  where  a consistent 

phase  jump  at  ky  = 0.6  was  evident  in  their  figures  for  9~  , Os'  , 

rll  12 

and  . Hence,  it  is  felt  that  the  structure  of  the  wave- induced 

r22 

flow  for  the  no-critical-layer  case  is  different  from  that  with  a 
critical  layer. 

The  data  in  general  agree  with  those  obtained  by  Chao,  et  al.  (1976), 
The  present  data  is  more  complete  than  theirs  because  the  wave-following 
measurement  in  this  study  covers  the  whole  boundary  layer  except  in  the 
immediate  vicinity  of  the  interface. 


5.2b  Critical  Layer  Flow  (Low  Critical  Layer,  U^/c  = 1.87  > 1) 

Figures  15  and  16  show  the  amplitude  and  phase  for  the  wave-induced 
velocities  u and  v , respectively,  for  U^/c  = 1.87  . As  shown  in 
Figure  15a,  the  amplitude  |u|  measured  in  the  fixed  frame  is  smaller 


than  that  measured  in  the  transformed  wave-following  frame.  This  is 
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opposite  to  what  we  observed  for  the  no  critical  layer  case.  This 

reversal  is  also  found  in  the  amplitude  distribution  of  |v|  shown  in 

Figure  16.  The  measurement  by  Hsu,  et  al.  (1977)  indicated  that  the 

phase  6^,  tends  to  be  90°  and  that  the  phase  9~  tends  to  be  270° 
u v 

near  the  interface.  Figures  15  and  16  clearly  show  this  tendency; 
however,  the  variation  in  9~  and  9~  at  this  higher  wind  speed  is  so 
rapid  that  our  measurement  does  not  reach  close  enough  to  the  interface 
to  reveal  accurately  the  phase  angles  at  the  interface. 

The  distributions  of  -u’v  are  shown  in  Figure  17.  They  have  a 
trend  similar  to  those  measured  by  Hsu,  et  al.  (1977)  and  by  Kendall 
(1970).  In  Hsu,  et  al.  (1977),  they  observed  a rapidly  increasing  region 
of  positive  values  of  -iiv  as  one  approaches  the  interface.  Such  a 
region  was  not  measured  here  because  it  is  too  close  to  the  interface 
at  this  wind  speed  for  our  probe  to  reach.  The  measurement  by  Kendall 
(1970)  indicated  that  this  characteristic  could  not  be  detected  with  a 
fixed  probe  even  at  lower  wind  speeds.  However,  a tendency  toward  a 
positive  value  in  - u v when  approaching  the  interface  is  seen  in 
Figure  17,  as  in  the  data  of  Kendall  (1970).  The  difference  in  the 

^ p u ru 

measured  value  of  -uv  (y)  and  - u v (y*)  is  the  value  of  f r)v 


is  positive  because 


9U 

9y* 


> 0.  Hence,  - u v (y)  > - uv(y*)  . This  feature 


is  apparent  in  Figure  17. 

Figures  18,  19  and  20  show  the  amplitudes  and  phases  for  r^  as 
functions  of  ky  at  this  wind  speed.  The  phase  jump  of  180°  at  the 
elevation  of  approximately  ky  = 0.A  is  clear  for  all  the  three  com- 
ponents of  9~  , 9~  and  8~  . This  is  consistent  with  the  measure- 

rll  r12  r22 

ment  by  Hsu,  et  al.  (1977).  Hence,  we  also  expect  that  the  amplitude  of 
r^j  should  have  a minimum  of  zero  at  approximately  ky  = 0.A  as  a 
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consequence  of  this  stress  reversal  there.  The  amplitudes  of  r^  are 
in  good  agreement  with  those  predicted  by  Hsu,  et  al.  (1978),  i.e., 


near  the  interface,  we  have  r^  w 2u'u'  kn  , u'u*  3n/3x  and 


~ 2uV  3n/9x  . Since  the  characteristics  of  r^  at  this  wind 
speed  are  similar  to  those  measured  by  Hsu,  et  al.  (1977)  for  U^/c  = 
1.57,  the  reader  is  referred  to  Hsu,  et  al.  (1977)  for  details  of  the 
structure  of  r^  . 


5.3  Probability  Distribution  a(y) 

Table  1 shows  the  result  of  the  joint  probability  between  the  u 
values  measured  by  the  two  probes  for  = 1.45  m/sec  . It  is  seen 
that  the  most  probable  occurrence  is  in  the  case  {u^  < c and 
u < c}  where  u is  the  horizontal  instantaneous  velocity  measured 

D SL 

by  the  old  probe  (higher  elevation)  and  u^  is  that  measured  by  the  new 

probe  (lower  elevation).  This  higher  occurrence  for  {u  < c and 

u,  < c)  is  expected  because  the  mean  velocity  u is  always  less  than  c . 
b 

However,  the  case  {u^  > c and  < c}  , which  is  favorable  for 
evaluating  the  momentum  and  energy  transfer  from  wind  to  waves  according 
to  the  modified  Miles'  theory,  happens  more  frequently  ttian  the  remaining 
two  cases.  For  < c , there  is  no  critical  layer  according  to 
Miles'  theory.  But,  near  the  free  stream,  turbulent  fluctuation  tends 
to  produce  an  instantaneous  velocity  profile  with  u greater  than  c 
at  the  outer  edge  of  the  boundary  layer  where  the  flow  is  strongly 
intermittent. 


1 
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2.92 


Table  2 shows  the  joint  probability  result  for  the  case  of  U » 

m/sec.  Now,  the  case  {u&  > c and  > c}  is  the  most  probable  case 

because  the  measured  u at  each  elevation  is  much  higher  than  c . 

Although  at  this  wind  speed  (U  > c)  is  different  from  the  previous  one 

(U  < c)  , the  favorable  case  of  {u  > c and  u,  < c}  is  still  higher 
00  a b 

than  the  remaining  two  cases.  If  we  include  the  turbulent  effect  and 
regard  the  statistical  result  of  this  time-varying  critical  layer  as 
a turbulently-mixed  critical  layer,  we  see  that  turbulence  tends  to  move 
the  mean  location  of  the  critical  layer  toward  the  location  where  the 
turbulence  has  its  highest  intensity. 

i 


i 
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6. 


CONCLUDING  REMARKS 


This  study  provides  a first  insight  into  the  difference  between 
the  wave-induced  flow  quantities  measured  in  fixed  and  wave-following- 
transformed  reference  frames.  The  two  measured  wind  speeds  cover  two 
different  flow  regimes  with  U < c and  U > c . The  wave-induced 

oo  co 

flow  for  these  two  flow  regimes  has  different  characteristics.  The 

phase  jump  for  r when  U > c may  be  crucial  in  the  wave  generation 
ij  00 

process  by  the  wind.  It  is  felt  that  the  two  wind  speeds  used  in  this 
study  were  not  very  efficient  in  defining  the  instantaneous  critical 
layer  relevant  to  the  modified  Miles’  theory.  For  a better  result,  a 
wind  speed  U just  above  the  wave  celerity  should  be  used. 
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TABLE  1 


Probability 

for  the  Occurrance  of  Critical  Layers  Above 
(Vjc  = 0.93) 

the  Interface 

A.  Wave-Following  Coordinate 

System 

y (cm)  | 

No.  of 
'ua  < c\ 

^ub  < c/ 

No.  of 
/ua  < c' 

\ub  1 ^ 

No.  of 

i 

No.  of 

Gil) 

a(y)  x 102 

2.6035 

184320 

0 

0 

0 

0 

3.1115 

184320 

0 

0 

0 

0 

3.8735 

184320 

0 

0 

0 

0 

4.8895 

184316 

0 

4 

0 

0.0022 

6.1595 

184310 

0 

10 

0 

0.0054 

7.4295 

184301 

0 

19 

0 

0.0103 

8.6995 

184306 

0 

14 

0 

0.0076 

10.6045 

184266 

0 

54 

0 

0.0293 

12.5095 

184206 

0 

114 

0 

0.0618 

15.0495 

183876 

0 

444 

0 

0.2409 

17.5895 

183072 

1 

1244 

3 

0.6771 

21.3995 

182020 

1 

2296 

3 

1.2478 

26.4795 

181727 

2 

2590 

1 

1.4068 

32.8295 

181431 

0 

2872 

17 

1.5674 

40.4495 

184257 

0 

63 

0 

0.0342 

B.  Fixed 

Coordinate 

System 

No.  of 

No . of 

No.  of 

No.  of 

y (cm) 

f"a  < c) 

(-  * ') 

(Ua  1 C) 

(■*-') 

a(y)  x io2 

\ub  < c/ 

\ub  > <7 

W < c/ 

\ub  lc/ 

2.6038 

— 

— 

— 

— 

— 

3.1115 

— 

— 

— 

— 

— 

3.8735 

— 

— 

— 

— 

— 

4.8895 

184320 

0 

0 

0 

0 

6.1595 

184320 

0 

0 

0 

0 

7.4295 

184312 

0 

8 

0 

0.0043 

8.6995 

184293 

0 

27 

0 

0.0146 

10.6045 

184229 

0 

91 

0 

0.0494 

12.5095 

184276 

0 

44 

0 

0.0239 

15.0495 

184159 

0 

161 

0 

0.0873 

17.5895 

183949 

2 

367 

2 

0.2013 

21.3995 

183286 

10 

1021 

3 

0.5610 

26.4795 

183981 

0 

332 

7 

0.1839 

32.8295 

183977 

0 

336 

7 

0.1861 

40.4495 

184319 

0 

1 

0 

0.0005 
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TABLE  2 


Probability  for  the  Occurrance  of  Critical  Layers  Above  the  Interface 

(Ujc  = 1.87) 


A.  Wave- 

-Following  Coordinate  System 

y (cm) 

No.  of 

r<C) 

\ub  < c/ 

No.  of 

( 

No.  of 

^b  < cj  \ 

No.  of 
ua  1 c\ 

ub  I c/ 

a(y)  x 102 

2.6035 

326 

499 

2009 

181486 

1.3607 

3.1115 

106 

381 

1017 

182816 

0.7585 

3.8735 

39 

144 

492 

183645 

0.3450 

4.8895 

5 

77 

88 

184150 

0.0895 

6.1595 

0 

0 

8 

184312 

0.0043 

7.4295 

0 

2 

6 

184312 

0.0043 

8.6995 

0 

0 

0 

184320 

0 

B.  Fixed  Coordinate  System 


y (cm) 

No.  of 

C; : :)' 

No.  of 
/u a < c 

W>c 

No.  of 

) C;::)  ( 

No.  of 
ua  1 c\ 

ub^C' 

ct(y)  x io2 

2.6035 

— 

— 

— 

— 

— 

3.1115 

— 

— 

— 

— 

— 

3.8735 

— 

— 

— 

— 

— 

4.8895 

6 

32 

113 

184169 

0.0787 

6.1595 

4 

1 

3 

184312 

0.0022 

7.4295 

0 

0 

0 

184320 

0 

a(y)  = Probability  of  [ua  < c , ub  c)  or  (ua  — c » ub  < 
Total  data  no.  = 184320. 

" — " means  that  it  cannot  be  measured. 
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Figure  15. 


Amplitude  and  Phase  Distributions  for  11  at  U^/c  » 1.87  . 
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Figure  1.  Schematic  of  Stanford  Wind,  Water-Wave  Facility. 
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Figure  4. 


Typical  Calibration  Curves  for  the  Hot-Film  Probes 
Al,  A2  — Films  1 and  2 of  Probe  No.  A (Old  Probe) 
Bl,  B2  — Films  1 and  2 of  Probe  No.  B (New  Probe) 
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Velocity  Profiles 


Figure  6.  Mean  Velocity  Profiles 


Figure  8. 


Profiles  for  Mean  Turbulent  Reynolds  Stresses  (U^/c  = 1.87)  . 
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Figure  9.  Amplitude  and  Phase  Distributions  for  u at  U /c  * 0.93  . 
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Figure  10 • Amplitude  and  Phase  Distributions  for  v at  U^/c 
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Figure  14.  Amplitude  and  Phase  E 


-5 


Figure  18.  Amplitude  and  Phase  Distributions  for  r at  U /c 
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